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ABSTRACT 

Slim disks have been received much attention because of the increasing evidence for the super-critical accreting 
objects. In this paper, we make an attempt to construct a unified model, in which the viscosity and the dimension- 
less accretion rate can span rather wide ranges. We replace blackbody radiation under diffusion approximation 
with a bridged formula, which accounts for both blackbody radiation and thermal bremsstrahlung in optically- 
thick and -thin cases, respectively. Thus this allows us to investigate the structures of and the emergent spectra 
from slim disks in a wider parameter space, covering transition regions from optically thick to optically thin. We 
show that there is a maximum transition radius, roughly /?tr.max/^g ~ 50 when M/Mc ^ 15. The emergent spectra 
from the unified model of the accretion disk have been calculated. A simple model of hot corona above the slim 
disk is taken into account for the hard X-ray spectrum in this paper based on Wang & Netzer (2003). We have 
applied the present model to the microquasar GRS 1915H-105, narrow line Seyfert 1 galaxies RE J1034H-396 and 
Akn 564. Our model can explain well the broadband X-ray spectra of narrow line Seyfert 1 galaxies, microquasars 
and possible ultra-luminous compact X-ray sources. The present model can be widely applied to the candidates 
of super-critical accreting objects. 

Subject headings: accretion, accretion disks — galaxies: active — X-rays: spectra 



1. INTRODUCTION 

It is believed that accretion disks around black holes are pow- 
ering in many kinds of celestial objects. The most important 
scale in the description of accretion disks is the critical accre- 
tion rate, defined by Mq = LEdi/rjc^ = 2.6 x IO'^Mbh/M© g/s, 
where the accretion efficiency 77 = 1/16 is for Schwarzschild 
black holes whose gravity effects are approximated by the 
pseudo-Newtonian potentials. The structures of accretion disks 
rely greatly on the mass of the black hole (Mbh), the accretion 
rate (M) and the viscosity (a), and hence the emergent spectrum 
from the disk. According to the dimensionless accretion rate 
m = M /Mc, the disk structures can be classified: (1) optically 
thin advection-dominated accretion flows (ADAFs; Narayan 
& Yi 1994), advection-dominated inflow-outflows (ADIOs; 
Blandford & Begelman 1999) and convection-dominated accre- 
tion flows (CDAFs; Narayan et al. 2000; Quataert & Gruzinov 
2000) when m < o?, (2) optically thick and geometrically thin 
disks (Shakura & Sunyaev 1973) if < m < 0.2, (3) slim disks 
(Muchotrzeb & Paczyriski 1982; Abramowicz et al. 1988; Chen 
& Taam 1993) once <m< 100. Chen et al. (1995) presented 
a unified model, but they assumed a Keplerian angular momen- 
tum distribution in the disk. The unified model to cover all these 
regimes has not been sufficiently understood, especially we still 
lack the emergent spectrum of such a model for the comparison 
with observations. 

There is growing evidence from current observations that the 
accretion rates of a significant fraction of celestial objects span 
quite wide ranges from very low to super-critical in X-ray bi- 
naries, microquasars and active galactic nuclei (AGNs). In our 
Galaxy, some ultra-luminous compact X-ray sources (ULXs), 
microquasars such as GRS 1915H-105 and SS 433, show large 
bolometric luminosities and high color temperatures (>1 keV) 



in their high states (Makishima et al. 2000; Watarai et al. 2001; 
Ebisawa et al. 2003). This provides strong evidence for the ap- 
pearance of slim disks in these objects. The Eddington ratios 
(L/LEdd) in AGNs and quasars have been reliably estimated for 
revealing the structures of the disks in AGNs and quasars, since 
the black hole masses can be estimated by the reverberation 
mapping relation (Kaspi et al. 2000) and the relation with their 
host galaxies (e.g. McLure & Dunlop 2001). It has been sug- 
gested that some of AGNs have close or super critical accretion 
rates (ColHn et al. 2002; Netzer 2003; Vestergaai'd 2003; Wang 
2003; Szuszkiewicz 2003), whose luminosities seemingly ex- 
ceed Lndd (if the estimations of black hole masses are reliable). 
Collin et al. (2002) have shown that about half of the 34 AGNs, 
of which the masses are determined by reverberation mapping, 
have super-critical accretion rates. Vestergaard (2003) extended 
the empirical relation of reverberation mapping in the sample 
of Kaspi et al. (2000) to the cases of high redshit quasars (see 
also McLure & Jarvis 2002). She found there are a signifi- 
cant number of quasars with super-Eddington luminosities, in 
either intermediate redshift range 1 .5 < z < 3.5 or high redshift 
range z > 3.5. Wiflott et al. (2003) found that the quasar SDSS 
Jl 148H-5251 at z = 6.41 is radiating at the Eddington luminosity 
with the estimated mass Mbh = 3 x lO^M©. Wang (2003) has 
also discovered a handful of the candidate super-critical accre- 
tors, using the limit relation between the black hole mass and 

FWHM(H^), Mbh = (2.9-12.6) x lO*" [ufwhm/1O^]*''^^M0. 
This relation is a result by combining the empirical reverbera- 
tion mapping and the features of the emergent spectra from slim 
disks. Although the spectra of 1 10 bright AGNs in the sample 
of Laor (1990) are fitted well by standard thin disks (Shakura 
& Sunyaev 1973), those super-Eddington objects can never be 
explained within the framework of the standard disk theory. A 
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wide range of dimensionless accretion rates should exist in var- 
ious kinds of celestial objects, but there is not a unified model 
of accretion disks to describe the varieties of accretion. 
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Fig. 1 . — Effective optical deptli (Teff) of tlie classical slim disk (see §2.1). 
The corresponding accretion rate is labeled. The effective optically thin region 
(Teff < 1) should be noted for 1 < m < 50. 

The significance of slim disks has been exhibited when they 
are applied to narrow line Seyfert 1 galaxies (NLSls; Wang 
et al. 1999; Mineshige et al. 2000; Wang & Netzer 2003; 
Kawaguchi 2003), Galactic black hole candidates (Watarai et 
al. 2000, Fukue 2000) and ULXs (Watarai et al. 2001) for the 
strong soft X-ray emission and the violent variability. However, 
slim disk models in above Uteratures suffer from the problem 
that the effective optical depth is smaller than unity in the inner 
disk region. The slim disks have been thought to work in these 
objects with high accretion rates, but the emergent spectra from 
the disks are not satisfactory enough. We list the models and 
spectra studied by far in Table 1, in which the main assumptions 
and features in these models are also given. The emergent spec- 
tra from the accretion disks of a wide range of rh remain open 
so far, especially only Wang & Netzer (2003) calculated the 
spectrum from the slim disk with hot corona based on the self- 
similar solution. We noted that most of the models assumed a 
low viscosity so that the transition to an optically-thin region 
does not appear in their models (see equation 1). Szuszkiewicz 
et al. (1996) calculated the emergent spectrum with the modi- 
fied blackbody radiation and explained the soft X-ray excesses 
claimed for some AGNs with slim disks. The vertical struc- 
tures and the continuum spectra have been computed by Wang 
et al. (1999) and Shimura & Manmoto (2003). Szuszkiewicz et 
al. (1996) presented the optically thin region with a small vis- 
cosity parameter a = 0.001 and the accretion rates 1 < m < 50. 
Kawaguchi (2003) also pointed out these regions in the accre- 
tion disks around super massive black holes (SMBHs). But they 
did not provide a reasonable treatment in the optically-thin hot 
region. In this paper, we show that a much wider region, where 
the effective optical depth Teff < 1, appears when the viscosity 
increases to q = 0. 1 for a stellar-mass black hole accretion (Fig. 
1). It can be seen that the transition radii from optically thick 
to thin locate at least several decades Rg (Rg = 2GM/c^, the 
Schwarzschild radius) for the accretion rates 10 < w < 50. In 
such large inner regions, the diffusion approximation for radia- 
tive transfer in the slim disks (Abramowicz et al. 1988) breaks 



down. We can not avoid the calculations of the emergent spec- 
tra from such disks having transition regions and hot corona, if 
connect the sUm disks with observations. 



Table 1 

A Brief Summary: Structures and Spectra of Slim Disk Models 
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no 


yes^ 


WZ99 


self-similar 


low 


no 


>1 


no 


no 


FOO 
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global 
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global 
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global 
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global 
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yes 
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self-similar 
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>1 
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global 
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4 
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Notes. — 

1: modified blackbody; 2: solving radiative transfer equation; 3: local 
blackbody approximation; 4: modified blackbody including Comp- 
tonization. 



Reference. — 

ASS: Abramowicz et al. (19SS); AOl: Artemova et al. (2001); B9S: Be- 
loborodov (199S); CT93: Chen & Taam (1993); C95: Chen et al. (1995); 
FOO: Fukue (2000); H96: Honma (1996); K03: Kawaguchi (2003); 
MOO: Mineshige et al. (2000); MPS2: Muchotrzeb & Paczyhski (19S2); 
SM03: Shimura & Manmoto (2003); SMA96: Szuszkiewicz, Malkan, & 
Abramowicz (1996); W99: Wang et al. (1999); WZ99: Wang & Zhou 
(1999); WFTOO: Watarai, Fukue, & Takeuchi (2000); WN03: Wang & 
Netzer (2003); ZTSOl: Zampieri, Turolla, & Szuszkiewicz (2001) 

This paper is organised by following: the basic formulations 
for the calculations of structures and spectra are given in §2, and 
the numerical results are presented in §3. We find that the tran- 
sonic location changes due to the presence of transition region. 
So does the structure in the inner disk region. We pay our atten- 
tions on how the emergent spectrum gradually changes with the 
accretion rate increasing from sub- to super-critical. The tran- 
sition to an optically-thin region modifies the emergent spec- 
trum, especially in the EUV and soft X-ray bands. We apply the 
present slim disk model with transition region to Galactic black 
hole candidates, one microquasar, and two NLSls sources in 
§4. The conclusions are laid in the final section. 

2. BASIC MODEL 

According to the model of Shakura & Sunyaev (1973), the 
effective optical depth will be less than unity at the radius 

RjRg^25a^/''m''/V/'\ (1) 

where we have neglected the factor (1 -r"'/^) (r = R/Rg and 
m = M/Mq). This formula shows that the transition covers 
the region where most of the gravitational energy is released. 
Within this region the radiation diffusion approximation does 
not work. Additionally, the transonic point is also located in 
this region. The transition radius /?tr is very sensitive to the 
accretion rate and the viscosity. We will focus on the global 
structure and the emergent spectrum from the slim disk with 
transition region. The equations are described as below. 
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2.1. Equations of slim disks with transition regions 

The basic equations of slim disks are taken from Muchotrzeb 
& Paczyriski (1982) and Abramowicz at al. (1988), but we deal 
with the radiative cooling more carefully. The vertical static 
equilibrium is assumed, namely H = B2Cs/i^K, where H is the 
half-height of the disk at the radius R, = \/P/p the sound 
speed, = \/GM]W /{l-R^/R) the Keplerian angular veloc- 
ity, P the total pressure and p the density. Unless stated, all 
physical quantities are referred to their equatorial plane values. 
We use the Shakura-Sunyaev viscosity prescription, = -aP, 
where is the viscous torque tensor and a is a phenomeno- 
logical viscosity parameter. The pseudo-Newtonian potential 
(Paczynski & Witta 1980) is adopted in our calculations for 
Schwartzschild black holes. 

The four equations of mass, radial momentum, angular mo- 
mentum and energy conservations control the structure of the 
sUm disk. The momentum equation in the radial direction is 



1^-(Q^- 
pdR 



ni)R+v^^=0, 



(2) 



where ur is the radial drift velocity of the accreting gas, fl the 
angular velocity. The specific angular momentum and the Kep- 
lerian angular momentum are defined by / and Ik = ^kR'^- The 
angular momentum conservation reads 

M(l-lin) = 4nR^HaP, (3) 

where /in is the eigenvalue of angular momentum at the inner 
disk edge. The mass accretion rate, M, is given by 

M = B52ttRT,-vr, (4) 

where S = 2pH is the surface density. 

In the classical slim disk model of Abramowicz et al. (1988), 
the viscosity parameter a is very low. This ensures the validity 
of the radiation diffusion approximation for the radiative cool- 
ing. However, the low viscosity only covers a very small param- 
eter space. When the viscosity increases, the density decreases 
so that the effective optical depth will be less than unity. In 
such a case the radiation diffusion approximation breaks down. 
In the regime of slim disk, we adopt a bridged formula 

4aT'^ (l-e-^-'OATeff 



Grad = ^3 ; 



(5) 



3(Tes + Tabs) C-^^f -|- (1 - e-^=«)Areff 

for the vertical radiation transport (see also Wandel & Liang 
1991; Chen et al. 1995). Here a is the Stefan-Boltzmann con- 
stant, T the mid-plane temperature, and the effective optical 
depth Teff = \/(''es + ''ahs)'''abs> whcrc wc take the electron scat- 
tering depth Tes = 0.34/?//, the absorption depth Tabs = 6.4 x 
lO^^p^T-^-^H for stellar -mass black holes, and Tabs = 1 -92 x 
10^^ p^T~^-^H for SMBHs due to bound-free processes (Czerny 
& Elvis 1987). This bridged formula works for a transition 
from optically thick to optically thin. The Comptonization am- 
plification factor A is adopted from Wandel & Liang (1991) 
for the simplicity. Compared with another type of bridged for- 
mula (Szuszkiewicz & Miller 1998; Narayan & Yi 1995), this 
bridged formula and the definition of Teff are consistent with 
those in spectral calculations described in below subsection. 
Heating due to viscosity, Q 

visj is cooled by local radiation, Grad? 
and global radial advection, Qadv, so the energy equation can be 
symbolically written as 

2vis = 2rad + 2adv, (6) 



where Qadv 



D M P c c 



-(12- l0.5l3)dlnT/dlnR+i4- 



ratio of gas to total pressure. The radial energy transport is ne- 
glected. The state equation is 



_ , 2rad(Tes + Tabs) 

pmn c 



(7) 



where /i = 0.617 is the mean molecular weight, kg the Boltz- 
mann constant, inn the mass of hydrogen. In this paper, we 
adoptfii = 0.67, /?2 = B3 = /^4 = 6 andfis = 0.5, just as Paczyriski 
& Bisnovatyi-Kogan (1981); Muchotrzeb & Paczyriski (1982); 
Muchotrzeb (1983). 

Two-temperature plasma is used in Wandel & Liang (1991), 
we assume the accreting gas is single-temperature over the en- 
tire disk, i.e. the energy transfer between the electrons and 
protons is very efficient. This assumption can be justified by 
comparing the timescales of free-free cooling and of Coulomb 
interaction between electrons and protons. The free-free cool- 
ing timescale is % = 2 x 10^p~[qT^-^ seconds and the timescale 
of the Coulomb interaction between electrons and protons is 
?Coui ~ 2 X 10~'^ p'^qT^-^ seconds, where pio = p/10""' and /g = 
T/IO^. When fcoui > fff, i-e. p < lO"'*" g/cm^ the plasma has 
two temperatures. We find that all the calculated p in this paper 
are consistent with the single-temperature assumption. 

2.2. Formulations for emergent spectra 

We employ the method described by Czerny & Elvis (1987) 
to calculate the emergent spectra from accretion disks, consid- 
ering the effect of opacity due to elastic electron scatterings 
and the effect of energy exchange in nonelastic scatterings (i.e. 
Comptonization). The photons are scattered before they escape 
from the surface of the disk (Rybicki & Lightman 1979), the 
emergent spectra are thus modified. However, Laor & Netzer 
(1989) pointed out that Comptonization is not important for 
disks around 10^ Mq black holes and L/Lem < 0.3. In slim 
disks, the transition region appears, where the free-free absorp- 
tion is not important but the scattering becomes essential. Wang 
& Netzer (2003) have shown this process is very important in 
slim disks, as confirmed by Kawaguchi (2003), see also §3.3 in 
this paper. 

The local radiation intensity 



IAr)=BUT,)f^(n,r), 



(8) 



where Bi, is the Planck function, Ts the temperature on disk sur- 
face, and fi, describing the departure from blackbody radiation 
including Comptonization is 

/.(7;,0 = /ff(l-/th)+C. (9) 

Here and the fraction of thermalized photons /th are given as 

2(i_e-2^«fO 



and 



ffi(Ts,r)-- 



f<h(Ts,r) = exp 



l + -\/l + Tes/Tabs,!.' 

ln(/vz.-7;/7»/) 
T3,ln[l + S,+ 16S2] 



(10) 



(11) 



3f})d}np/d'hiR is the dimensionless advection factor, /? is the 



respectively, where the electron's dimensionless temperature 
Ee = AksT^meC^, and t^* = (Tes + Tabs,i/)/(l + Teff). /ff rcstores 
the modified blackbody radiation for Teff ^ 1 , and describes the 
thermal bremsstrahlung from an optically thin medium for low 
Teff. The Comptonization normalization constant is calculated 
by 

2>kBTs f^-Klv/hvdv 



C = 



(12) 
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The surface temperature Ts is reached by the balance between 
the radiated energy and the dissipated gravitational energy 

/■OO 

erad= / nUR,%)dv (13) 

This simple method for allows us to avoid the calculations of 
the vertical structure of the slim disk, and the accuracy of the 
spectrum is enough for the present aims. The total emergent 
spectrum is derived by integrating /^(r) over the radius, 

L^ = 2 I TTl^iR)2TTRdR. (14) 

The hard X-ray emission from the hot corona above the slim 
disk is also calculated based on the detailed description in Wang 
& Netzer (2003). The most important characterized parameter 
/, the fraction of the gravitational energy released in the hot 
corona, is assumed in our calculation. However, this factor / is 
confined by the fact that the hot corona should keep, at least, as 
optically thin medium, otherwise it will be rapidly cooled. This 
limit can be roughly estimated by Tes < 1 at the last stable ra- 
dius with the assumption that the radial velocity is light speed. 
Wang & Netzer (2003) showed 

(til \ ^ 
To) • ^^^^ 

This limit shows that the hot corona becomes weak with the 
increase of the accretion rate. Indeed the hard X-ray spectra 
from microquasars and NLSls show such a property. In our 
calculations, we thus neglect the feedback of hot corona to the 
slim disk. This is justified because the factor / is small for a 
high accretion rate. However, a detailed treatment including 
this feedback for the case of intermediate accretion rate will be 
carried out in a separate paper. 

3. NUMERICAL RESULTS 

The equations (2-7) can be reduced to two first-order differ- 
ential equations about two variables, p and T. We use the shoot- 
ing method to solve these equations with appropriate boundary 
conditions. We assume an outer boundary at lO'' Rg, where 
Shakura-Sunyaev solutions (the gas pressure-dominated and 
Kabs ^ Kes) are available. At this boundary, we firstly choose 
the derivatives of p and T obtained from Shakura-Sunyaev disk 
solutions, substitute them into the two differential equations, 
then calculate the values of p and T by iterations. Although 
they are very close to Shakura-Sunyaev values, the new val- 
ues of p and T allows the inward integration to run on the 
right way. The general relativistic effects are simulated by the 
pseudo-Newtonian potential near the horizon of black holes. 
Fourth-order Runge-Kutta method is employed to integrate the 
two differential equations to the inner disk edge R\n, where Rin 
is taken as a free parameter as other authors did (e.g. Mineshige 
et al. 2000; Wang & Netzer 2003). We use the torque-free con- 
dition at the inner boundary. The shooting method determines 
the eigen angular momentum /„,, provided the solution satisfies 
the regularity condition at the critical point close to the tran- 
sonic radius. 

After we obtain the global transonic solutions, the basic 
physical quantities in the accretion flow and its emergent spec- 
trum are given for stellar black holes (IOMq) and SMBHs 
(IO^Mq) in this section. We show in detail how the transition 
to an optically thin region appears with the change of accretion 
rate. 



3.1. Transonic solutions 
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Fig. 2. — Angular momentum (/) and total pressure (P; top panel), ratio of 
radial velocity to sound speed (vr/cs; bottom panel) for m = 1 (solid lines) 
and m = 10 (dashed lines). The total pressiu-e of ni = 1 is shifted upward by 
one logarithmic unit. The transonic locations are marked by solid points. The 
grey lines are the results of the classical slim disks calculated in this paper (the 
same for Figs. 3-5, 9, 10 and 12). 
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Fig. 3. — Location of the transonic point (Ss) in function of the accretion 
rate (m) and the viscosity parameter (a; labeled for each line). 

The detailed transonic situations are analysed here, and they 
are compared with the cases of slim disks without the correct 




treatment of transition regions. Fig. 2 presents the exact dis- 
tributions of the angular momentum (/), totoal pressure (P) and 
radial velocity (wr) in the vicinity of transonic position, with the 
viscosity parameter a = 0.1 and m = 10. The cases of m = 1 and 
10 are both exhibited. From this figure, it can be seen that the 
curves of m = 1 of our model are almost the same as those of the 
classical slim disk model. Comparing with the classical slim 
disk, in the slim disk with transition regions, the angular mo- 
mentums shift upward, total pressures decrease and the ratios 
of radial velocity to sound speed increase at the innermost ra- 
dius, due to the effects of optically thin region. In the inner disk 
region where the effective optical depth is smaller than unity, 
the diffusion approximation for blackbody radiation is invalid, 
and the radiative coohng becomes inefficient. Thus the temper- 
ature is higher than that of the classical slim disk, where local 
blackbody radiation is assumed in the entire disk. At the same 
time, the lower total pressure resulting from the decreased radi- 
ation pressure due to Teff < 1, leads to the smaller viscous tensor 
(Try oc P). Therefore, the curve of angular momentum is flatter 
(see equation 3) and the flow becomes more dense. The accre- 
tion flow reaches the sonic point earlier than the classical slim 
disk, since the sound speed = (P / pf -^ is much smaller and 
the ratio wr/ Cs becomes larger. This explains why the transonic 
positions shift outwards, compared with those in Abramowicz 
etal. (1988). 

As well known, the accretion flow has four types of transonic 
cases, according to which kind of pressure (gas or radiation) is 
dominant at the sonic point and the magnitude of the viscos- 
ity parameter. For low a, it is total pressure that pushes the 
flow pass through the sonic point. While for large a, the action 
of viscosity drives the flow across the transonic point. Fig. 3 
shows the dependency of the transonic radius on the accretion 
rate and the viscosity parameter for a black hole of 10 Mq. It 
is apparent that the transition to an optically thin region moves 
the transonic point out, when the accretion rate is in the range 
1 < m < 100 and the viscosity parameter is a > 0.01. In other 
parameter space comprising rh and a, the transonic location in 
our result is almost the same as that of the classical slim disk. 

The transonic flows with transition regions of high viscosity 
can be compared with the flows of low viscosity. First, they 
are both radiation pressure dominated for super-critical accre- 
tion rates. Second, they are different in some properties. One 
of them is manifested in the optical depth. The accretion flow is 
almost entirely optically thick for a small viscosity parameter, 
e.g. a = 0.0001. While the transition region appears for a large 



a (if a is large enough, the inner flow can be purely optically 
thin and gas-pressure dominated). Another distinct feature lies 
in the position of the sonic point. For low a, the flow cannot 
acquire enough velocity to transonic until it is near the marginal 
bound orbit (2 Rg), since the sound speed is relatively high and 
the viscosity-driven process is not efficient. On the contrary, the 
sound speed of the flow with large a is small, the sonic point 
locates far outside of the last stable orbit (3 Rg). Therefore, we 
conclude that the transonic situation of the accretion flow with 
transition region is, to some extent, between the case of low 
viscosity and that of high viscosity in the classical sUm disk. 

The S-shaped curves for different parameters are shown in 
Fig. 4. It is found that they are strongly affected by the transi- 
tion region when the viscosity is high. For the case of stellar- 
mass black holes, the upper ADAF and the middle unstable 
branches change much. The upper branch due to advection is 
modified, and the unstable branch is compressed. The appear- 
ance of the transition to an optically thin region lowers the radi- 
ation efficiency, hence the advection becomes more important. 
In the accretion flows of SMBHs, the S-shaped curves are less 
altered, as shown in panels (c) and (d) in Fig. 4. It is easy to 
understand the curves are not seriously influenced in the outer 
region, since there the optical depth is large and it tends to be 
the classical sUm disk. The distorted S-shaped curves may have 
observational features in stellar black hole accretions. The tran- 
sition timescale from the upper to the lower branches, ti , will 
be shorter than that of the classical slim disk. For example with 
the same parameters, the difference at /? = 5/?g has a factor of 
2-4. But the transition timescale from the lower branch to the 
upper, f|, is unchanged and greater than f|. This would be of 
interest in the explanation of X-ray transient objects. The vari- 
ation of state-transition timescales in stellar black holes is also 
one of the crucial characteristics of the transition to optically 
thin regions in sUm disks. 

3.2. Global disk structure 

Fig. 5 provides the physical quantities such as the half- 
thickness, density, midplane temperature (including surface 
temperature), effective optical depth and the ratio of gas pres- 
sure to total pressure for slim disks of both stellar-mass black 
holes and SMBHs. The results of the classical slim disks are 
also shown in the figure with grey lines. The transition region 
influences greatly the structure in inner disk region, especiaUy 
for the accretion flow around a stellar-mass black hole. 

First is the half-thickness. For aU the super-critical accretion 
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Fig. 5. — m-dependency of the aspect ratio (HjR), density (p), midplane temperature (7) and surface temperature (Ti), effective optical deptli (r^fi), ratio of gas 
pressure to total pressure (/3) for black holes of (a) m = 10 and (b) m — 10* with a = 0. 1. The solid, dashed, dot-dashed, dotted and dot-dot-dot-dashed lines are for 
?n = 1, 3, 10, 20 and 50, respectively. The surface temperatures % are shifted down by 0.5 logarithmic units for clarity. 



rates 1 < m < 50, < 1, indicating the flow is geometrically 
slim. The thickness increases with the accretion rate, which 
is consistent with the result of the classical slim disk. At the 
outer boundary radius, the disk is always geometrically thin, 
and the solutions can be approximated by Shakura-Sunyaev so- 
lutions. Because ofH = B2Cs/flK, the lower sound speed results 
in the reduction of half-thickness, in comparison with the clas- 
sical slim disk model. The reduced height favors the trapped 
photons to escape from the disk. This shape of the accretion 
flow is somewhat analogical with the furmels in thick disks. 



With the emergence of the transition region, the temperature 
and density of the disk are both larger than those in the classical 
slim disks. There are still considerable radiation fluxes in the 
innermost region after the flow passes through the sonic point, 
which is pointed out by Watarai et al. (2000). The effective 
optical depth has also been shown in Fig. 5. It should be noted 
that there is a large zone of Teff < 1 when m > 3, e.g. more 
than 30 R„ (15 R^} for the disk of a black hole of 10 (10*") solar 
masses. The lowest effective optical depth reaches 0.03. More 
important is that, in the accretion flows of both stellar-mass and 
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supermassive black holes, there are wider regions of at least 
100 Rg showing Teff ~ 1 for a moderate super-critical accretion 
rate 10 < m < 50. In such large transition regions, blackbody 
radiation is not adequate to describe the radiative transfer. Ac- 
cordingly, a new radiation formula valid for any optical depth, 
as adopted in the present work, is necessary. 

When the accretion rate tends to very large, for example 
m> 10, in the inner disk region the advected energy dominates 
over the surface cooling, i.e. gadv ^ Qrad, the flow will have 
self-similar behaviors (Wang & Zhou 1999). We find that the 
global solution indeed shows the self-similar signatures: the ra- 
dial velocity ur oc total pressure P oc r~^/-^, the density 
p oc r"^/-^, the temperature T oc and Teff oc r^/'* weakly 
depends on the radius, except in the vicinities of the inner disk 
edge. Watarai & Mineshige (2003) have also shown these prop- 
erties. 
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Fig. 6. — The dependence of the transition radius on the accretion rate for a 
fixed viscosity a = 0.1 in stellar black holes. The black points are the results 
from our slim disk model. There is a maximum radius of the transition region. 
The solid line is equation (1) for Shakura & Sunyaev (1973) disk. The dotted 
is the fitting line by the least square method. 

Fig. 6 shows the dependence of the transition region on the 
accretion rate. We find the relation can be well fitted by an 
analytical formula 

Rtr/Rg = 95 3Sa°''^rfP-^^exp (-0.1m° ^) . (16) 

There is a maximum transition radius Rtr.max/Rg ~ 50 at wi 15. 
This is clearly different from that in the standard disk model, 
i.e. Rtr oc m^/^^ (see equation 1). Our calculations show that 
there is no transition radius when rh < 0.5, and it drops quickly 
if m > 30. The maximum radius of the transition region can 
be understood from the S-shaped curves in Fig. 4. Near the 
turning point between the upper and the middle branch, the sur- 
face density is the lowest and at the same time the very high 
temperature makes the absorption opacity be smaU. The father 
from this area, the larger effective optical depth due to the in- 
crease of either the surface density or the absorption opacity, 
or of both. This is an interesting result since the intermediate 
super-critical accretion has the largest transition radius, where 
the gravitational energy is mostly released. 

We briefly discuss the viscosity effect on the disk struc- 
ture. Fig. 7 is the a-dependency of the global disk structure. 
In the slim disks, a has only a Uttle effect on the disk half- 
thickness. In our model large a such as 0.3, reduces much 



height, increases the density and temperature in the inner re- 
gion, which differs with that discussed by Kawaguchi (2003). 
At the outer radius, the situation is similar as the result in 
Kawaguchi (2003): the larger viscosity, the lower density, mid- 
plane temperature and also the optical depth, but the higher sur- 
face temperature. This is because of more efficient angular mo- 
mentum transfer outwards and more trapped advection energy 
in the flow. 
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Fig. 7. — a-dependency of the physical quantities for m = 10 and m = 10. 
Solid, dashed, dot-dashed and dotted lines are for a = 0.3, 0.1, 0.01 and 0.001, 
respectively. The surface temperatures (Ts; grey lines) are shifted downward 
by 0.2 logarithmic units. 

Fig. 8 is the dependence of the transition radius on the 
viscosity a. It can be found that the results from our calcu- 
lations are not in agreement with the Shakura-Sunyaev disk, 
even at a small viscosity {a < 0.01). This is because that the 
Shakura-Sunyaev disk breaks down when m > 0.2. The de- 
pendence of the transition radius on a becomes very strong, 
roughly Rtr oc a" "^' for a > 0.01 and oc a"-22 for a < 0.01, 
while Rtr the Shakura-Sunyaev disk. The high vis- 

cosity drives the accreting gas to fall into the black hole more 
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rapidly, leading to a low density region. This region then be- 
comes absorption weak, enlarging the transition region. Since 
the radial motion and the advection effect is neglected in the 
standard disk, its prediction does not hold when either the ac- 
cretion rate or the viscosity is large. 

In the end, we should pay attention that the properties of the 
effective optically thin region discussed here are completely 
different from opticaUy-thin ADAFs. Firstly, the flow is only 
effective optically thin. The volumn density and the surface 
density are still larger than those of optically-thin ADAFs. The 
Thomson scattering is very important and the Comptonization 
is significant. That the optical depth becomes effectively thin is 
because of the transparancy to photon absorption at a high tem- 
perature. In optically-thin ADAFs, the plasma density is very 
low and the Thomson scattering optical depth is rather small. 
Secondly, the temperature in the present slim disk model is not 
so high as that of optically-thin ADAF (10^"' K). Thirdly, flie 
pressure is still radiation-pressure dominated. While for the 
ADAF, it has merely gas pressure, the radiation pressure can be 
neglected. However, we believe that slim disks (optically-thick 
ADAFs) are possible to transit to optically-thin ADAFs in inner 
region, once the viscosity is large enough. We leave this topic 
in the later paper. 
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Fig. 8. — The dependence of the transition radius on the viscosity for a 
fixed accretion rate m = 10 in stellar black holes. The black points are the re- 
sults from our slim disk model. The solid line represents the equation (1) from 
Shakura & Sunyaev (1973) model. The dotted line is the fitting Une. 

3.3. Emergent spectra 

We have shown that the global structures of the slim disks 
with transition regions are quite different from those slim disks 
in Abramowicz et al. (1988). With the presence of transition 
region, only the theoretically predicted spectrum based on the 
global structure is of great significance to the observations. 

The emergent spectra of slim disks have been extensively 
studied. Szuszkiewicz et al. (1996) were the first to calculate 
the spectrum with modified blackbody radation, and explained 
the soft-X excesses of quasars. Wang et al. (1999) calculated 
the vertical structures and emergent spectra from slim disks. 
Two prominent features appear: the saturated total luminosity 
and the constant cutoff energy. Mineshige et al. (2000) and 
Watarai et al. (2001) used purely blackbody radiation approx- 



imation, though it is not adequate in inner disk region. Wang 
& Netzer (2003) and Kawaguchi (2003) employed the method 
of Czerny & Elvis (1987) to consider Thompson electron scat- 
terings and Comptonization in accretion disks. However, since 
they did not treat the phenomena of optically thin, their results 
are only available for low viscosity parameter. As for the case 
of high viscosity parameter, the local radiation in spectral calcu- 
lations tends to thermal bremsstrahlung when the optical depth 
is small, while local blackbody radiation that holds only for 
optically thick medium is assumed to determine the disk struc- 
ture. This inconsistent treatment between structure and spectral 
calculations leads to the problem of temperature invertion (Cz- 
erny & Elvis 1987), namely, the calculated surface temperature 
is even higher than the midplane temperature in the accretion 
flow (see the grey lines in the middle panels of Fig. 5; In the 
classical shm disk, is much larger than T at several /?g). The 
reason is that the blackbody approximation predicts too much 
energy loss than the realistic radiation in optically-thin flows. 
Thus, the radiated flux is heavily overestimated at high ener- 
gies. 

We have solved the problem in this paper. By adopting the 
bridged formula which can describe the radiation energy for 
both optically-thick and -thin cases, the calculated surface tem- 
perature is not larger than the midplane temperature, as shown 
in Figs. 5 and 7. Except that the Thompson electron scatter- 
ing modifies the radiation to deviate blackbody, the inelastic 
scattering also exchange the energy between photons and elec- 
trons. The Comptonization effect has been stressed in Wang 
& Netzer (2003) and Kawaguchi (2003). We can give a more 
detailed estimation of the importance of this process. The ef- 
fective Compton parameter 

3;, = S,Max(r4,T,'2), (17) 

where = Tes/Max(reff, 1), the absorption opacity is taken as 
Kabs,!. = 1 .4 X W-^ pT-'^-^x-\\-Q'') {x = hv/ksT), and 30 times 
value for SMBHs. Fig. 9 shows the parameter for different 
photon energies, accretion rates and black hole masses. When 
> 1, Comptonization can not be neglected and the emergent 
spectrum will be strongly modified. Fig. 9 clearly shows how 
important the Comptonization process is for disks with different 
parameters. We can draw a conclusion that the Comptonization 
effect is very important in the inner region of super-critical ac- 
cretion flows. Saturated Comptonization process is inevitable 
for less massive black hole disk (lO^^M©). This conflrms 
the conclusion of Wang & Netzer (2003). The Comptoniza- 
tion is less important for a 1O^M0 black hole disk unless the 
disk has highly super-critical accretion rate, which agrees with 
Laor & Netzer (1989). The higher accretion rate, the more 
important Comptonization process. We adopt the scheme de- 
scribed in §2.2 to calculate the spectra from slim disks, taking 
into account the Comptonization effect by a simple but practical 
method. 

Figs. lOand 11 show the details of the emergent spectra from 
the sUm disks with transition regions. We would hke to stress 
that the present calculations have two advantages: 1) we include 
the transition region; 2) the model covers wide ranges of the ac- 
cretion rate and the viscosity. Thus the present model involves 
all the known disk models (except for the optically thin ADIOs 
and CDAFs). Fig. 10 clearly shows the spectra from the ac- 
cretion disks of which the central black hole masses are 10 and 
10^ Mq. We start the accretion rate from sub-critical (m = 0.3) 
to super-critical (m = 50). We find the spectrum shifts toward 
high energy with the increasing accretion rate. This is due to 
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(a) log R/R^ (m=10) 



(b) log R/R^ (m=10'') 



(c) log R/Rg (m=10'') 



Fig. 9. — Effective Compton parameter y* in slim disks with transition regions. In each panel, the solid lines correspond tox = hhi/ksT = 0.2, 0.4, 0.6, 0.8 and 

1 .0, respectively, from bottom to top. 



the increases of temperature and of the importance of Comp- 
tonization with the accretion rate. The spectra show very strong 
humps at high energies for both stellar and supermassive black 
holes. The spectral flux is proportional to the accretion rate for 
the sub-critical accretion flow. The distinct spectral character- 
istic of super-critical disks is its broadness and flatness. The 
cutoff energies of the spectra almost keep constant when the 
accretion rate m is above 10, and the total radiated luminosity 
from the disk tends to be saturated. The constant cutoff energy 
was first found in Wang et al. (1999), which is caused by ad- 
vection. The feature of the saturated luminosity has also been 
found in Abramowicz et al. (1988), which is caused by the ad- 
vection effect: most of the dissipated energy is swallowed by 
the central black hole. The maximum cutoff energy is about 
100 keV (4 keV) for a black hole of 10 (10^) solar masses. The 



saturated luminosity is about several Eddington limit (^ 10^^ 
and 10"*"* erg/s for 10 and 10* Mq black holes, respectively). 
This extreme slim disk with very high accretion rate can be ex- 
plained by the behaviors of self-similar solution (Wang & Zhou 
1999; Wang & Netzer 2003). The features of the broad and 
flat spectrum, constant cutoff energy and saturated luminosity 
are very useful in explaining the X-ray spectra of microquasars, 
ULXs and narrow line Seyfert 1 galaxies. 

For a comparison with the spectra from the classical slim 
disks, we draw them as grey lines in Fig. 10. We use the 
same computation scheme for the spectra from the classical 
slim disks. We find that the emergent spectra are quite dif- 
ferent from those based on the present model with transition 
region, when the accretion rate is above the critical value. The 
cutoff energy in the classical slim disk is unrealisticaUy higher 
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Fig. 10. — The emergent .spectra from .slim disks with transition regions for (a) m = 10 and (b) m = 10*. From bottom to top, the accretion rates are 0.3, 1, 3, 10, 
20 and 50, respectively. The grey lines are the results for the classical slim disks, of which the calculated surface temperature is higher than the midplane value in 
the optically thin region. Thus the radiation fluxes at high energies are overestimated. 





Fig. 1 1 . — a-dependency of the emergent spectra from slim disks with transition regions. The solid, dashed, dot-dashed and dotted lines are for a = 0.001, 0.01, 
0.1 and 0.3, respectively. The accretion rates are labeled in each panel. 



than the present model. Thus the present model provides an 
improved version for the slim disk. 

We investigate the effect of the viscosity on the spectrum. 
For all the accretion rates, larger viscosity parameter causes 
the spectrum to extend to higher frequency. This is the con- 
tribution of the more significant elastic electron scatterings and 
Comptonization. From Czerny & Elvis (1987), the critical tem- 
perature is oc m^/'^OT'/"Q"^/", above which the modifica- 
tion due to electron scattering is important. If a increases, the 
transition region will enlarge and the surface temperature will 
increase. In such a flow the Comptonization effect becomes 
stronger. Thus the emergent spectrum is broader and flatter 
with the higher viscosity parameter. Usually the cutoff energy 
strongly relies on both the accretion rate and the viscosity pa- 
rameter. But in a saturated accretion system, it is only affected 
by a. This may be utilized to determine the unknown viscos- 
ity parameter through X-ray observations of ultra- luminous ob- 
jects. 



Finally, with the detailed calculations, we verify the impor- 
tance of photon trapping in the super-critical flows. If the dif- 
fusion timescale 

. g(Tes + Tabs) 

'diff= J (io) 

C 

is longer than the accretion timescale 

4cc = — (19) 

in a flow, the photons can NOT escape from the disk at the same 
radius. This causes the so-called photon trapping. Such a case 
has been recognized in early time (Katz 1977; Begelman 1978), 
but the quantitative features have been discussed by Wang & 
Zhou (1999); Ohsuga et al. (2002, 2003) and Shimura & Man- 
moto (2003). 

Let us compare the two timescales. Fig. 12 shows face ^ ?dyn 
in the entire r < 10^ region for any rh < 50. This means that 
the viscosity works for the disk, transporting the angular mo- 
mentum and dissipating the gravitational energy. When m > 1, 
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in the inner region fdiff > face, the photons are hardly radiated 
out so that trapped in the accretion flow. In this trapped region, 
energy transport in the form of advection is the main cooling 
mechanism. In this sense, the slim disk is often called op- 
tically thick ADAF. After comparisons, we conclude that the 
photon trapping effect is still significant, but it is not so se- 
vere as that in the classical slim disk. This is mainly because 
of the greatly reduced disk height. The present model does 
not treat the photon trapping carefully, but see Ohsuga et al. 
(2002, 2003), who treated the photon trapping with the given 
disk structure. Wang & Zhou (1999) have shown that the ra- 
diated luminosity is of 4 x 1O-'^(Mbh/1OM0) erg/s for an ex- 
treme super-critical disk, which is as small as one third of the 
Eddington luminosity. In fact, in physics the photon trapping 
effect makes it difficult to estimate the real accretion rate from 
the saturated observed spectrum, when the black hole has a 
super-critical accretion rate. One of the goals in this paper is 
to investigate the physics of the transition region. For a mod- 
erate super-critical rate m < 10, the photon trapping effect is 
not so strong, our present model does work for our applications 
in next section. The model coupling the photon trapping ef- 
fects is just in preparation. Considering the photon trapping ef- 
fect, we believe that the super-critical accretors in microquasars 
and NLSls are more than the number we know at the present. 
Kawaguchi (2004) reaUzed the role of the self-gravitation in 
the optical spectrum of NLSls. We did not consider the self- 
graviation in the present work. The general relativistic effects 
are not included either Future work will solve the equations 
of the general relativistic disk with transition region, as well as 
calculate the spectrum from the disk. 




log R/Rg 

Fig. 12. — Accretion timescale (face; ttiick lines), diffusion timescale {tag; 
thin lines) and dynamic timescale ((dyn). The soUd, dashed, dot-dashed and 
dotted Hnes are for m = 1, 10, 20 and 50, respectively. The grey lines are the 
results of the classical slim disks. 

4. APPLICATIONS 

The most prominent characters of the spectra from slim disks 
are the strong humps and the constant cutoff frequency. When 
L/LEdd > 0.3, the standard accretion disk model breaks down 
(Laor & Netzer 1989) and the slim disk is available. Our present 
model covers a rather wide range of accretion rates. The dis- 
tinguished soft X-ray humps in some objects are explained by 



the self-consistent emergent spectra based on the corresponding 
global disk structures. In this section, we apply our model of 
slim disks to the microquasar GRS 1915H-105 and two NLSls: 
RE J1034-H396 and Akn 564. 

4.1. Microquasars: GRS 1915+105 

Some of galactic black hole candidates, especially micro- 
quasars and ULXs, are regarded to be super-critical accretors 
(Watarai et al. 2000, 2001). The broad band X-ray spec- 
tra of microquasars and ULXs have been observed extensively 
with RXTE and ASCA (Makishima et al. 2000; Mizuno et al. 
2001; Ueda et al. 2002). The observed spectra are fitted well 
with multicolor blackbody spectra in the framework of stan- 
dard optically-thick accretion disk. However, the multicolor 
disk (MCD) spectrum has some shortcomings. The spectral 
hardening factor due to Comptonization in the disk atmosphere 
is assumed to be a constant (usually 1.7; Shimura & Maimioto 
2003) over the entire disk. Our calculations show Comptoniza- 
tion is only crucial in inner region. Additionally, MCD is based 
on the standard disk theory, the viscous dissipation is thought 
to be radiated efficiently in the local. This endows an upper 
luminosity limit (L < 0.3Lndd; Laor & Netzer 1989) to its use 
scope. Obviously, the luminosities of ULXs and other bright 
objects are beyond this limit. Therefore, the fitting parameters 
of MCD model have difficulties to connect with the physics. 
First, the X-ray estimated black hole mass with MCD is smaller 
than the precise dynamic measurement by optical observations. 
Second, the fitted color temperature is too hot to be produced 
by the standard disk, i.e. the problem of too hot disk (e.g. King 
& Puchnarewicz 2002; Ebisawa et al. 2003). Lastly, tiie MCD 
model cannot account for spectral variations. The five periods 
of spectral fitting of IC 342 source 1 (one ULX) show that the 
black hole mass varies by a factor of 2 (Ebisawa et al. 2003). 
This is unreasonable. 

Ebisawa et al. (2003) found the spectra of Galactic super- 
luminal jet sources can be explained with standard Kerr disk 
model, while ULXs need slim disk model. Sobolewska & Zy- 
cki (2003) stated that the relativistic effects of Kerr model can 
only partially account for the apparently complex soft X-ray 
spectra of GRS 1915H-105. In this subsection, we show that the 
emergent spectrum from our slim disk model, considering the 
effects of deviated blackbody radiation, Comptonization and 
the transition to optically thin, is quite adequate to describe the 
soft spectral component of GRS 1915-H105 when it is in high 
state. 

The measurements of the black hole mass and the distance 
in GRS 1915-1-105 are quite reliable, so it is able to apply 
the present model to this object. The black hole mass Mbh = 
(14 ± 4)Mq is deduced from the direct measurement of the or- 
bital period and mass function of GRS 1915H-105 (Greiner et 
al. 2001). The inclination is i = 70°, the distance T> = 12.5 
kpc (Mirabel & Rodriguez 1994). We adopt the recent values 
i = 66° and I? =1 1 kpc (Fender et al. 1999). 

GRS 1915-1-105 stays almost permanently in a soft spectral 
state, although it has strong variability of both flux and energy 
spectta. RXTE data have been extensively studied by Belloni et 
al. (2000). The lightcurves are divided into 12 classes accord- 
ing to the color-color diagrams and the count rates. Each class 
can be decomposed into only 3 basic spectral states: A, B and 
C. Sobolewska & Zycki (2003) attempted the modeling of these 
three states, using the data of PC A and HEXTE onboard RXTE. 
In the framework of MCD, they found that the X-ray spectra of 
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A and B (corresponding to the relatively higher luminosity) are 
complex: only one Comptonized disk plus another model such 
as Comptonized disk, blackbody or Kerr disk can give an ac- 
ceptable fitting. They derived the accretion rate M = 7.5 x lO'^ 
g/s, namely m ^ 0.3 for the spectral state B. We choose the 
same observation data to do fitting with the present model. For 
a comparison, we employ the absorption cross section given by 
Morrison & McCanunon (1983) to correct the spectrum from 
our slim disk. The effective hydrogen column density is taken 
as 6.2x10^2 (Sobolewska & Zycki 2003). The spectral 
energy distribution is shown in Fig. 13. The dotted Une is 
the spectrum from the slim disk described in this paper. The 
Kerr MCD spectral component of Sobolewska & Zycki (2003) 
is also shown by the grey dashed line. The hard tail is modeled 
by the thermal Comptonized spectrum (dashed line) from the 
corona above the sUm disk. With our sUm disk model, the best 
fitted parameters are a = 0.01, rh = 2.1 and / = 0.03. It can be 
seen that the theoretical spectrum coincides with the observa- 
tion very well in the whole X-ray band. With these parameters, 
the transition radius is ~ 3.4 7?g. In this case, the difference 
from the classical slim disk is not distinguishable. But it is clear 
that our model is more self-consistent in the inner disk region. 
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Fig . 13. — X-ray spectra of GRS 1 9 1 5+ 1 05 after an absorption by the hydro- 
gen column density of 6.2x lO^^cm"'^. Solid line is the fitting spectrum with 
our transition slim disk (dotted line) plus the corona (dashed Une). The RXTE 
observation data and the multi-color Kerr metric disk spectrum (grey dashed 
line) are taken from Sobolewska & Zycki (2003). 

Since our slim disk covers rather wide ranges of accretion 
rate and viscosity parameter, it can be also appUed to other stel- 
lar black holes, such as galactic black hole binaries, superlumi- 
nal jet sources (microquasars) and ULXs. We have also suc- 
cessfully fitted the X-ray spectra of GRO J1655-40 (Tomsick 
et al. 1999), LMC X-1 (Gierlinski et al. 2001). Their accre- 
tion rates are ^ 0.5 Mq- Ebisawa et al. (2003) applied the slim 
disk spectra to IC 342 source I (see also Waterai et al. 2001 for 
other ULXs). Since the basic quantities of ULXs, e.g. the mass 
and the distance, arc unknown, the spectral fitting is a little dif- 
ficult. Once we implement the present model into XSPEC, it 
will be possible to model many other sources, especially XTE 
J 1550-56 and IC 342 source 1 suspected having super-critical 
accretion rates. This will help us understand more physics in 
ULXs. Making use of the wide validity of the present model, 
we can understand the process of state transition, of which the 
transition to an optically thin region exists. 



4.2. Narrow line Seyfert 1 galaxies 

Narrow line Seyfert 1 galaxies (NLSls), as a special sub- 
group of Seyfert 1 galaxies, are characterized by their unusual 
properties: the H/3 line is narrow (FWHM H/3 < 2000 km/s) 
and relatively weak (less than a third of the intensity of the 
[OIII]A5007 line), and the optical Fell fines are very strong. 
It is revealed that the excess and the tendency of the spectrum 
to flatten at low energy, have been found in several NLSls, 
e.g. PKS 0558-504 (O'Brien et al. 2001), Mrk 766 (Boiler et 
al. 2001), PG 1244-H026 (Fiore et al. 1998), PG 1404H-226, 
PG 1440H-356, PG 1211H-143 (George et al. 2000), IRAS 
13224-3809 (Vaughan et al. 1999), IH 0707-495 (Boiler et 
al. 2002, Dewangan et al. 2002), RX J1702.5H-3247 (Gliozzi 
et al. 2001) and Ton S180 (Turner et al. 2001a, 2002). Bep- 
poSAXmAASCA observations of Akn 564 and RE J1034H-396 
clearly show these extreme properties including a soft X-ray 
hump with a very flat spectrum (vFi, = const.; Comastri et al. 
2001 ; Puchnarewicz et al. 2001 ; Turner et al. 2001b; VignaU et 
al. 2003; Romano et al. 2004). Other unusual X-ray properties 
are discussed in Brandt et al. (1994); Comastri et al. (2001); 
Pounds, Done, & Osborne (1995) and Collinge et al. (2001). 
Earlier suggestions attributing the soft excess to blends of emis- 
sion fines are not supported by recent A^CA, XMM-Newton and 
Chandra observations (e.g. Turner, George, & Netzer 1999; 
Puchnarewicz et al. 2001; Turner et al. 2001a, 2002b; CoUinge 
et al. 2001). Chandra shows that there is not any feature of 
emission fines in soft X-ray band (Collier et al. 2001). The 
model of less massive black hole with high accretion rate is 
preferred for its very prominent soft X-ray hump (Boiler et al. 
1998; Wang & Netzer 2003; Kawaguchi 2003). In fact the fea- 
ture of soft X-ray hump vF^ = const is a natural explanation of 
super-critical accretion disk (Wang & Netzer 2003). With the 
black hole mass estimated by the reverberation relation, Wang 
& Netzer (2003) found half of NLSls have L/Leaa > 1 in the 
sample of Vef on-Cetty (200 1 ) . We befieve that NLS 1 s are good 
candidates for the applications of the present model. 

As an appUcation of the present model, we make an attempt 
to explain the observed soft X-ray humps of Akn 564 and RE 
J1034H-396. RE J1034H-396 shows a very clear appearance of 
flattening of soft X-ray spectrum from logv = 16.5 to 17.2 in 
the uL^ plot (Comastri et al. 2001; Puchnarewicz et al. 2001; 
Collinge et al. 2001; Turner et al. 2001). Here we suggest 
that this feature reflects the power of the slim disk in the two 
objects. 

RE J1034+396: The ultra-soft X-ray radiation of RE 
J1034-H396 was first discovered by Puchnarewicz et al. (1995). 
The spectrum shown in Fig. 14 shows a clear soft X-ray 
hump extending from \ogv = 16.5 (the lowest observed fre- 
quency) to 17.2. The data were taken from Puchnarewicz et al. 
(2001). They explained the strong soft X-ray spectrum of RE 
J1034H-396 based on the standard accretion disk model (Czemy 
& Elvis 1987). The resulting accretion rate is m ^ 0.3-0.7, 
which beyonds the scope of the standard disk model. They also 
suggested a high incUnation j'obs = 60-75° in order to lower the 
luminosity. 

The upper panel of Fig. 14 shows the comparison of our 
model with the observation data. The FWHM of H/3=2000 
km/s and the luminosity L5100 = 10"*"^ erg/s at 5100 A, the mass 
of the black hole is estimated as Mbh = 2.25 x lO^M© by the 
empirical reverberation mapping relation (Kaspi et al. 2000). 
We assume the inclination angle i = 0°, since no accurate esti- 



Slim Disk and Application to Micro-quasars and Narrow Line Seyfert 1 Galaxies 



13 



mation has been made in the literature. The viscosity a = 0.2 is 
required from the fitting. We obtain the accretion rate m= 1.2. 
The factor / = 0.07 is found, which is within the limit of equa- 
tion (15). This result is in a rough agreement with that of self- 
similar solution in Wang & Netzer (2003). 




log v (Hz) 

Fig. 14.— Broad-band spectra of REJ 1034+396 (top) and Akn 564 (bot- 
tom). The X-ray spectra are both observated with BeppoSAX, and the data are 
taken from Puchnarewicz et al. (2001) and Comastri et al. (2001). The dotted 
lines are the slim disk spectra. The hard tail is modeled by the Comptonized 
spectra (dashed lines) from the corona above the slim disk, using our self- 
consistent disk/corona model. The transition radii to the optically thin regions 
are 11.5 and 15.5 Rg for RE J1034+396 and Akn 564, respectively 

Arakelian 564: Pounds et al. (2001) estimated that the 
black hole mass in Akn 564 is ^ 10^ Mq and the accretion 
rate m « 0.2- 1 from the break frequency of the power spec- 
trum of the light curves. The multiwavelength spectra of Akn 
564 have been extensively studied by Romano et al. (2004). 
They found that Mbh = 4 x lO^M© and m s=; 1. Wang & Net- 
zer (2003) provided a simple method to estimate the black 
hole mass if the accretion rate is super-critical, Mbh = 2.8 x 
10^ (i^L,y/lO'*'*erg/s)M0, which yields the black hole mass 
2.0 X 10*" Mq (Romano et al. 2004) for Akn 564. The current 
data has a few factors of uncertainty of Mbh- The inclination 
angle i = 27.4° (Ballantyne et al. 2001), which was derived 
from the fitting ASCA spectra with the ionized reflection disk 
models of Ross & Fabian (1993). The observed spectral energy 
distribution of Akn 564 is taken from Comastri et al. (2001). 
The comparison of the present model with the observations pro- 
vides Mbh = lO^M©, m = 1.1, a = 0.32 and / = 0.021. The hot 
corona is relatively weaker than that in RE J1034-I-396. 

The fitting results indicate that the two NLSls are moderate 
super-critical accreting sources. In this case, the transition re- 
gion plays great role in the emergent spectrum. The transition 



radii, where Teff = 1, are - 1 1 .5 and l5.5Rg for RE J1034H-396 
and Akn 564, respectively. This is just the main X-ray emit- 
ting zone. Therefore, although the transition region influences 
little S-shaped curves (Fig. 4), it has important effect on the 
spectrum from the disk. Our slim disk model provides more 
accurate treatment about this transition region, so it has a wider 
application field than the previous slim disk model. As stated by 
Pounds & Reeves (2002), the soft X-ray humps are stiU poorly 
understood since we do not know what drives the shape of soft 
X-ray spectrum. We need more observations. The future work 
will focus on the spectral fitting in a sample so that we can 
discover the primary parameters controlling the shape of soft 
X-ray spectrum. 

It is noted that the mass accretion rates for the three objects, 
GRS 1915H-105, RE J1034-H396 and Akn 564, are of 1 - 2 crit- 
ical accretion rates. For this range of accretion rate, the shm 
disk with a viscous law have a large thermally unstable region. 
These unstable regions are responsible for the observed rapid 
and giant variabilities from microquasars and NLSls. To de- 
scribe the structures accurately, the time-dependent slim disk 
model with transition regions is needed to explore the variabil- 
ities of these super-critical accretors. We leave this topic in a 
future research. 

5. CONCLUSIONS 

The detailed structure and spectral calculations of the slim 
disk with transition region (the effective optical depth Teff < 1) 
have been presented. We have shown that there is a quite 
large transition region in the slim disk. This region be- 
comes larger with the increase of viscosity parameter. We get 
an empirical formula for the transition radius, i.e. Rtr/Rg = 
95.38a'^^'m''^^exp (-O.lm'*'). The maximum transition ra- 
dius reaches ^ 50 Rg at m ~ 15 for a fixed viscosity a. In 
the slim disk regime, the transition radius is Ra/Rg oc a"^' for 
a fixed accretion rate (e.g. m = 10) and a > 0.01. The transition 
region affects not only the global disk structure, but also the 
emergent spectrum. The transonic location moves outward and 
the local radiated spectrum departs from blackbody radiation, 
because of the presence of the transition region. In compari- 
son with the classical slim disk, we have made more accurate 
treatment with the optically thin region, in the calculations of 
both structure and spectrum. The present model covers a rather 
wide parameter space. We have also shown how the soft X- 
ray spectrum changes from sub-critical, moderate super-critical 
to the extreme super-critical accretion rate. We calculate the 
Comptonization parameter for the slim disk, and show that the 
Comptonization is very essential for stellar-mass black holes 
and less massive black holes (< 10*"^ Mq), whereas it is not 
very important for more massive black holes with the masses 
larger than lO' Mq. This shows that the Comptonization pro- 
cess play a key role in the formation of spectra in luminous 
objects. The present model can be widely applied to micro- 
quasars, ultra-luminous X-ray sources and narrow line Seyfert 
1 galaxies. 
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